On the black hole from merging binary neutron stars: how fast can it spin? 



Wolfgang Kastaun/ Filippo Galeazzi,^'^ Daniela Alic/ Luciano Rezzolla/ and Jose A. Font^ 

^ Max-Planck-Institut fur Gravitationsphysik, Albert-Einstein-Institut, Potsdam, 14476, Germany 
^Departamento de Astronomiay Astro flsica, Universitat de Valencia, Dr. Moliner 50, 46100, Burj as sot (Valencia), Spain 

The merger of two neutron stars will in general lead to the formation of a torus surrounding a black hole whose 
rotational energy can be tapped to potentially power a short gamma-ray burst. We have studied the merger of 
equal-mass binaries with spins aligned with the orbital angular momentum to determine the maximum spin 
the black hole can reach. Our initial data consists of irrotational binaries to which we add various amounts of 
rotation to increase the total angular momentum. Although the initial data violates the constraint equations, 
the use of the constraint-damping CCZ4 formulation yields evolutions with violations smaller than those with 
irrotational initial data and standard formulations. Interestingly, we find that a limit of J/M^ 0.89 exists for 
the dimensionless spin and that any additional angular momentum given to the binary ends up in the torus rather 
than in the black hole, thus providing another nontrivial example supporting the cosmic censorship hypothesis. 
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Introduction. A long-standing open question in general rel- 
ativity is how close the angular momentum J^^ of a black hole 
(BH) created in astrophysical scenarios can get to the the- 
oretical maximum, which for an isolated BH is , with 

' ADM' 

M^j^j^ the ADM mass, and under very general conditions is 
A/(87r), with A the horizon area Q. So far, no violation 
of the cosmic censorship hypothesis, i.e., no formation of a 
naked singularity, has been found in numerical simulations 
evolving generic initial data satisfying the dominant energy 
condition (but see also |2|). This question has been studied 
for the collapse of supermassive stars, e.g., |3|, neutron stars 
(NSs) with J^s/^ns ^ ^1^1' merger of binary NSs 

(BNSs), e.g., IMLl of BH-NS binaries, e.g., |8 |, and of BH 
binaries, e.g., In all cases investigated, the final BH spin 
is below the critical value and semi-analytical estimates seem 
to indicate this is true for all configurations ifTOlfTTIl . 

Determining the final spin produced in the merger of BNSs 
is not a mere academic question as the rapid rotation of the 
BH is a key ingredient in all of the models in which the BNS 
merger is thought to lead to a jet formation and a short gamma- 
ray burst |[T2l - [T4B . Because the energetics of the emission will 
depend sensitively on the BH spin, an accurate measure of 
the maximum value attainable can help set upper limits on 
the efficiency of the process (see 1 15 1 for a recent discussion 
and a list of references). Despite its importance, this issue 
has not been addressed yet and the main reason is the lack 
of constraint- satisfying initial data for spinning NSs. Early 
attempts to construct such initial data, e.g., |[I6|[T21, did not 
have a satisfactory solution of the Euler equation according 
to 1, 1 8 ] . Only recently, initial data for BNSs with spins has 
been computed |[T9l , but no evolutions have been reported yet. 

Here we follow a novel approach, which consists of setting 
up constraint- violating initial data by adding a rotational ve- 
locity field to self-consistent irrotational initial data, and then 
evolving with a recently developed constraint-damping CCZ4 
formalism [20]. We show that this choice reduces the con- 
straint violations quickly to a level which is even smaller than 
the one encountered in BNS evolutions of irrotational initial 
data with standard formulations. The drawback of this simple 
method is that the artificial spin-up also introduces oscillations 
of the star and eccentricity of the orbit. This would be unac- 



ceptable for modelling gravitational-wave emission, but it is 
adequate to assess the influence of the NS spin on the BH. 

Binary Initial Data with Spin. The best way to decide what 
are realistic expectations for the spin in merging BNSs is to 
look at the available observations. Although the distribution 
of NS spins is still unknown, NSs at birth are expected to have 
spins in the range 10-140 ms ll2T1l . and dimensionless spin pa- 
rameters J^^/M^^ ~ 1 should be attainable only by millisec- 
ond pulsars that have been spun-up by accretion |22|. It is 
however unlikely that the old NSs of a binary system which is 
about to merge have periods less than 1 ms, and all the obser- 
vational evidence supports this conclusion. We know, in fact, 
that the period of the fastest known pulsar in a BNS system, 
J0737-3039A, is 22.70 ms, yielding a dimensionless spin pa- 
rameter J^^/M^^ ~ 0.05 |23|. In view of these considerations 
we restrict our analysis to binaries — 0.2 < J^^/M^^ ^0-3, 
thus well within the range of realistic possibilities. 

To construct initial data for a binary system of spinning 
NSs, we modify the velocity field of the irrotational equal 
mass solutions computed using the LORENE code |24|. We 
set = {1 — s) wl -\- sft X X, where w'^ = /u^ denotes 
the coordinate velocity of the fluid in the star, the fluid 
4-velocity, the orbital angular velocity vector, is the 
original irrotational velocity field, and s is a free parameter 
we tune to increase/decrease the spin of the star J^^ (s = 
corresponds to a purely irrotational binary, 5 = 1 to a corotat- 
ing one). Once a spin-up/down is introduced, we parametrize 
a binary in terms of the additional ADM angular momentum 
it has with respect to the irrotational model, AJ^j^j^. Note 
we compute J^^^ using Eq. (68) of 1241 because this form is 
more robust against constraint violations. Two remarks should 
be made on our prescription for w. First, it corresponds to 
stars with spins aligned with the orbital angular momentum, 
which is the most relevant case for our goals, since it leads 
to the maximum increase of total angular momentum. Sec- 
ond, it guarantees that the residual vector field is orthogonal 
to the density gradient, and hence the flow is adapted to the 
deformed shape of the stars in the binary. However, because 
the hydrostatic equations are no longer satisfied with the new 
velocity field, it introduces oscillations in the stars as they are 
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evolved. The relative variation of the central rest-mass den- 
sity stays below 1% for the irrotational models, is < 8% for 
s = 0.85, and < 15% for the extreme case of 5 = 1.2. 

Using the above method, we construct sequences of in- 
creasing spin, each with a fixed baryon mass, from 1.625M0 
to I.9OIM0 for each star, and the same initial separation of 
45 km, which corresponds to roughly 3 orbits until merger for 
the irrotational case (up to 4 for the fastest spinning cases). 
The stars obey initially a polytropic equation of state (EOS) 
P{p) = Kp^ with r = 2 and = 123.6 in units in which 
G = c = Mq = 1. During the evolution, however, we use 
an ideal-fluid (F-law) EOS |5 1 with F = 2. Since we are in- 
terested in the maximum spin of the BH, we focus on binaries 
that will collapse promptly as any intervening long-lived hy- 
permassive neutron star (HMNS) would just extract additional 
angular momentum from the system. Hence, we consider only 
models with total baryon masses well above the maximum one 
for a nonrotating star, which for our EOS is 2.0 Mq. We in- 
crease the spin frequency up to 55% of the Kepler limit, which 
corresponds to 5 = 1.2, but larger values are possible in prin- 
ciple. As in our previous work (e.g., |25 1), we use adaptive 
mesh-refinement techniques 1 26 1 with 6 levels during inspiral, 
two of which follow the stars, and a 7th finest level activated at 
the time of collapse. The outer boundary is located at 756 km, 
where Sommerfeld radiation boundary conditions are applied. 

CCZ4 and BSSNOK formulations. We have already men- 
tioned the importance of using a formulation of the Einstein 
equations that damps rapidly the violations of the constraints 
caused by our modifications of the initial data. To achieve 
this, we employ the recently proposed conformal and covari- 
ant Z4 formulation (CCZ4) 1 20 1 and couple it to the equations 
of relativistic hydrodynamics. We have carried out a system- 
atic investigation of this new formulation in comparison with 
the more standard BSSNOK one |27|, including the evolution 
of a stable isolated NS, the collapse of an unstable star to a 
BH, and the merger of BNSs in quasi-circular and eccentric 
orbits. For lack of space these tests will be presented in de- 
tail in a forthcoming longer publication and here we restrict 
the discussion to the constraint-damping properties of CCZ4. 
The comparison reveals not only the advantages of CCZ4 over 
BSSNOK. It also provides a way to estimate how much con- 
straint violations are tolerable before the influence on the dy- 
namics becomes significant, by looking at the difference be- 
tween two simulations of the same system with very different 
constraint violations. More specifically, we can start assessing 
the role played by the initial constraint violations with the ex- 
treme case of an eccentric binary as created from quasicircular 
initial data in which we reduce artificially the linear momenta 
by 15%. Doing this induces very large constraint violations 
in both formulations, although the L2-norm of the Hamilto- 
nian constraint in CCZ4 evolutions always becomes smaller 
by a factor 5-10. Not surprisingly, the dynamics of the sys- 
tem shows large differences in the phase evolution and in the 
trajectories, demonstrating that the initial constraint violations 
are too high to yield meaningful results. On the other hand, for 
quasicircular binaries, the orbital trajectories for CCZ4 and 
BSSNOK agree quite well, with a phase error of 4%. At the 
time of merger, the mass and the spin of the BH agree within 



0.6% and 0.4%, which is less than our estimate of 1% for the 
numerical errors. Again, we find that the use of CCZ4 reduces 
the L2-norm of the Hamiltonian constraint, on average by one 
order of magnitude. In particular, a very sharp decrease is 
seen during the first millisecond (cf. left panel of Fig. [T]). Un- 
der similar but not identical conditions, Ref. 1 28 1 has reported 
a larger decrease for an alternative conformal formulation of 
the Z4 system, Z4c. Without a direct comparison it is difficult 
to assess why it results in a smaller violation. It may also be 
due to the improved treatment of the outer boundaries in 1281 . 

A similar behaviour is also observed for NSs that are spun- 
up/down. This is shown in the left part of Fig.[T} which reports 
the evolution of the L2-norm for a number of binaries with 
different degrees of spin-up/down, when evolved with either 
the CCZ4 or the BSSNOK formulation. The initial violation 
is obviously larger for larger changes of spin, but in all CCZ4 
evolutions the violation is reduced by about an order of mag- 
nitude after 1 ms, becoming less than the one from BSSNOK 
evolutions with constraint-satisfying initial data. Afterwards, 
the constraints exhibit the usual behavior resulting from the 
interplay between constraint-violating numerical errors and 
the action of the damping terms in the CCZ4 system |20|, re- 
sulting in a series of maxima and minima. This behaviour 
should be contrasted with the one of the BSSNOK formu- 
lation, where the constraints instead remain large during the 
whole evolution. For the latter, the errors tend to accumulate 
at each point in space, causing the stars to leave behind lumps 
of constraint violations. This is shown in the right part of 
Fig. [T] which displays the Hamiltonian constraint violation at 
different stages of the evolution for the case s = — 1, both for 
CCZ4 and BSSNOK. The comparison highUghts that while 
the constraint violations are damped/propagated away in the 
CCZ4 evolution, they remain static in the BSSNOK one. 

Interestingly, even with the most severe spin change of the 
s = —1 binary, the violation with the CCZ4 formulation 
is smaller by a factor ~ 5 than the one obtained with the 
BSSNOK formulation for irrotational binaries (red solid line), 
which represents the best-case for that formulation. There- 
fore, since in all CCZ4 evolutions the constraint violations fall 
below the irrotational case with BSSNOK, we conclude that 
after the first millisecond the errors due to constraint viola- 
tions are comparable with the differences between BSSNOK 
and CCZ4 in the irrotational case. Furthermore, to assess also 
the influence of the large violation during the first millisec- 
ond we note that for the case s = —1, the difference in BH 
mass, spin, and merger time between BSSNOK and CCZ4 
evolutions are 0.7%, 3%, and 10%, respectively. Since CCZ4 
evolutions always have smaller violations after the first mil- 
lisecond, we conclude that the error in the BH mass and spin 
due to constraint violations to be less than 1%, and the total 
error, including the numerical one, to be less than 2%. 
Final BH spin. We turn now to discuss our physical results 
and their as trophy sical implications. We first note that for 
the lightest binary (M^^ = 1.63 Mq), only the irrotational 
model forms a BH promptly, while a HMNS is formed for 
s = 0.5, 1.2. This is not unexpected, but highlights that if the 
total mass is close to the critical mass for a prompt collapse, 
the spin of the NSs can have a strong impact on the dynam- 
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FIG. 1. Left. Time evolution of the L2-norm of the Hamiltonian constraint H over the whole domain, excluding the AH, obtained using CCZ4 
and BSSNOK formulations. The symbols mark the AH formation. Right. Snapshots in the orbital plane for s = — 1, at times (a) t = 0, (b) 
t — 1.48 ms, during inspiral, (c) t — 3.48 ms, shortly before horizon formation, (d) t — 4.46 ms, shortly after. The left half of each panel 
shows the L2-norm of H as obtained with CCZ4 evolutions, the right half with BSSNOK. The red solid lines mark the trajectories of the NSs. 
In all cases the binaries have M^^ = 1.78 Mq. 



ics of the merger and thus on the gravitational- wave signal. 
For the heavier binaries, instead, we observe the immediate 
formation of a BH and a torus. Quasi-stationarity is reached 
soon after merger, with the BH mass and spin changing less 
than 0.4% during the last 50 M^^ for all simulations, and with 
negligible mass accretion. These conditions are essential to 
use the isolated-horizon formalism to measure the properties 
of the BH. 

To measure the contribution to the total angular momentum 
of the fluid orbiting outside the BH we monitor the quantity 
Jp = - /^^ d^Vn^^^T^^, where d^V, n^.T^u denote, re- 
spectively, the portion of the time slice outside the apparent 
horizon (AH), the proper volume element, the unit normal to 
St, the energy-momentum tensor, and = 0^ is the basis 
vector of the spherical coordinate system obtained from the 
Cartesian one used in the simulations. Jp reduces to the an- 
gular momentum in the Newtonian limit, but is not conserved 
unless is a Killing vector (in which case it is conserved even 
for non-axisymmetric flows 1291 ) and coincides with the Ko- 
mar angular momentum for stationary axisymmetric space- 
times. Once the BH becomes stationary, Jp is the angular mo- 
mentum of the matter in the torus. 

On the other hand, to measure the angular momentum of 
the BH, we use the isolated-horizon formalism 1301 [3T1 , com- 
puting the integral J^^{t) = (Stt)"^ J^d^V^^R'' K^j, on the 
AH surface At, where <l>^, i?^, K^j^ are, respectively, an (ap- 
proximate) axial Killing vector on the AH, the unit normal to 
At on and the extrinsic curvature. Surprisingly, the sum of 
Jp and Jgj^ at the end of the simulations agrees with Jp at the 
time of AH formation to better than 2.2% on average, with a 
variance ^ 0.5 % during the whole time from the formation of 
the horizon until stationarity. This behaviour can be explained 
using Eq. (45) in 1321 . which indicates that the time derivative 
of this sum is small if (j)^ becomes approximately a Killing 
vector of the spacetime and the gravitational-wave emission 



from the AH is small. Both conditions are satisfied at the time 
of horizon formation to a degree sufficient for us to use Jp at 
this time to estimate the final total angular momentum. 

The results of our simulations are summarized in the left 
panel of Fig.[2j which reports with filled symbols the dimen- 
sionless spin of the BH, J^^/M^^. Note that M^^^ is the hori- 
zon mass 1 31 1 computed from J^^ and the horizon area using 
the Kerr formula, so that the extremal case SttJ^^^ = | 1 | is 
equivalent to J^^^ = M^^ . The data is presented as a function 
of the (dimensionless) additional initial angular momentum 
per star, ^AJ^^^/M^^, relative to the corresponding irrota- 
tional model. Values of AJ^^^/M^^ > correspond to an 
increase of total angular momentum, while AJ^^^/M^^ < 
to NSs with antialigned spins. For comparison, we also plot 
the estimate for the dimensionless spin of PSR J0737-3039A. 

As it is natural to expect, the BH spin increases with the 
initial total angular momentum of the binary, and for a binary 
like PSR J0737-3039A it is not much larger than for an irro- 
tational one. binary. Surprisingly, however, the growth is not 
constant and actually saturates for AJ^^^/M^^ 0.2, quite 
independently of the initial (baryon) mass in the binary. The 
largest BH spin we obtain is J^^/M^^ = 0.888 ± 0.018 for 
Mj^3 = 1.78 Mq and s = 1.2. At the same time, we also find 
that the combined angular momentum of the BH and of the 
torus, Jtot = Jbh + Jj^ continues to increase for all values of 
s used. More precisely, for s < 1, the increase is essentially 
linear, with Jtot = J^°^ + V^Jadm^ where r] = 0.55 (0.60) 
for = 1.78 (1.9O)M0 and J^^^ is the value of Jtot for 
the irrotational binary. This means that roughly half of the ad- 
ditional angular momentum has been radiated as gravitational 
waves. For sufficiently large spin-ups, a significant fraction of 
the remaining angular momentum is transferred to the torus, 
which therefore acts as the channel absorbing the excess an- 
gular momentum and limits the spin-up of the BH. This is 
reported in Fig. |2] which shows that Jp (empty symbols) is 
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FIG. 2. Dimensionless spin of the BH J^^/Ml^ (filled symbols) and 
of the torus J^/M^^ (empty symbols) versus the additional initial an- 
gular momentum A J^j^j^. The error-bars include discretisation errors 
and the influence of the constraint violations. 



tion of the Einstein equations from generic initial conditions 
leads to a BH formation rather than a naked singularity. 

We will conclude with a few remarks. First, it is reasonable 
that the material in the torus will eventually be accreted onto 
the BH, transferring angular momentum and further increas- 
ing the BH spin. However, this will happen on dissipative 
timescales which are much longer and thus not relevant for 
the central engine of gamma-ray burst, which should be ig- 
nited promptly with the merger. Second, because magnetic 
fields are expected to further decrease the angular momentum 
in the system they are not particularly relevant in our consid- 
erations, which focus on the maximum spin in a prompt BH 
formation. Finally, although the BH spin seems to increase for 
irrotational unequal-masses 1331 . the spin could also decrease 
depending on the amount of mass and spin in the torus. In any 
case, we expect an upper limit will be found also for generic 
binaries and this will be the focus of our future work. 



very small for antialigned spins and irrotational binaries, but 
increases rapidly with the initial stellar spin. 

Of course it is difficult to provide a simple explanation for 
this nonlinear behaviour, although one might speculate that 
spinning up the NSs will provide matter with sufficient angu- 
lar momentum to then produce a torus. Overall, however, our 
results provide another nontrivial example that, as proposed 
by the cosmic censorship hypothesis, a self-consistent evolu- 
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